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, Abstract. Here, we consider damped two-components Bose- Einstein 

^ ' condensates with many-body interactions. We show that, when 

^ , the external trapping potential has a double-well shape and when 

' ' ' the nonlinear coupling factors are modulated in time, hysteresis 

^ : effects may appear under some circumstances. Such hysteresis 

' ' phenomena are a result of the joint contribution between the ap- 

i, pearance of saddle node bifurcations and damping effect. 

^ ■ PACS number(s): 05.45.-a, 02.30.Oz, 03.65. Ge, 03.75.Lm 

> ■ 

■ Atomic Bose-Einstein condensates (BECs) at zero temperature are 

^ , ! described by means of nonlinear Schrodinger equations of the type 

(1) ^n^ = Ho^lj + gM''^P + 9M'^P + ... 

I ■ at 

^ i where Hq represents the Hamiltonian of a single trapped atom and the 

^ I nonlinear term \^\'^^, r = 1,2, . . ., is the (r -|- l)-body contact poten- 

CnI ■ tial with coupling factor gr+i [I]. In fact, BECs strongly depend by 

^ . interatomic forces and the binary coupling term (yf2|V'PV' usually rep- 

^ I resents the dominant nonlinear term; when the higher nonlinear terms 

Q I are neglected then equation ([1]) takes the form of the well-known Gross- 

^ ■ Pitaevskii equation [2J. The coupling factor of the binary nonlinear 

>■ , term is given by (72 = N'2nh'^a/m where m is the mass of the atoms, 

I A/" is the total number of particles of the condensate, and a is the scat- 

^ ' tering length; for higher nonlinearity terms some expressions of the 

coupling factor have been recently proposed [3]. Recent experiments 
[1] have shown that the scattering length a can be changed, and, in 
fact, the 2-body coupling factor g2 can be tuned to be zero in the case 
of polar molecules in optical lattices driven by microwave fields [5]. In 
such a case, the 3-body and, in general, (r -|- l)-body interaction be- 
comes significant [6l [7] and thus the Gross-Pitaevskii equation becomes 
inadequate in order to describe BECs. 

The basic properties of BECs with many-nody interactions described 
by equation ([T]), where many nonlinear terms are simultaneously con- 
sidered, are far to be well understood. In order to understand some 



Date: June 17, 2010. 



1 



2 



ANDREA SACCHETTI 



fundamental features of BECs with many-body interactions in this pa- 
per we follow the approach proposed by [S], that is we are interested 
to the (r + l)-body intereaction in its own. Therefore, we restrict 
ourselves to the basic nonlinear Schrodinger equation representing the 
properties of exactly (r + l)-body contact interaction of BECs at zero 
temperature 



of a single atom in dimension n with trapping potential V{x), as well 
as on the (r + l)-body coupling factor Qr+i- 

It is worth mentioning also the fact that equation ([2]) with nonlinear- 
ity corresponding to the power-law \'ip\'^^ where the parameter r takes 
any positive real value, is used in other contexts, including semicon- 
ductors [9] and nonlinear optics [ini [HI [E]- Furthermore, even if in 
most of the applications the parameter r takes only integer and posi- 
tive values, here we take that r can assume non integer values too, as 
considered in [13]. 

Since the wavefunction if) is assumed to be normalized to 1 then the 
coupling factor gr+i depends, in addition to the physical parameters of 
the problem such as the scattering length and the mass of the parti- 
cles of the condensate, on the total number AT of the particles of the 
condensate. 

Recent experiments have shown that the total number M of particles 
that participate to the condensate can be adiabatically modulated in 
time by means of a suitable time-dependent combination of optical 
and magnetic forces [15] . Thus, we can consider the case where 
the nonlinear coupling factor in equation ([2]) is a given function which 
slowly depends on the time t 



Another way to produce a time-dependent coupling factor consists of 
tuning the scattering length [16]. In a previous theoretical papers 
by Pelinovsky, Kevrekidis and Frantzeskakis [T7] the Gross-Pitaevskii 
equation with a periodically varying nonlinearity coupling factor has 
been considered and it has been shown a good agreement between 
solutions of the averaged and full equations. 





which depends on the Hamiltonian 




dr+l '■— gr+l{t) ■ 
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In this paper we show that the modulation of the nonhnearity cou- 
phng factor may give rise to hysteresis phenomenon for two- components 
BECs, where the the external trapping potential V{x) has a double- 
well shape [H]. In fact, hysteresis effects have been already seen in 
rotating BECs, in particular the number of vortices appearing in a ro- 
tate BEC depends on the rotation history of the trap, in addition to 
the number of vortices initially present in the condensate [19] ; see also 
the theoretical analysis in |2U] . 

In particular, here we show that, in the semi classical limit of h small 
enough, for BECs equation ([2]) with power law l^pl"^^ in a double well 
trapping potential and with a slow modulation, with respect to the 
beating period between the two wells, of the nonlinear coupling factor 
gr+i{t) then an hysteresis effect appears provided that r is bigger than 
the critical value 



It is worth mentioning the fact that this result holds true for both 
attractive (e.g. gr+i < 0) and repulsive (e.g. gr+i > 0) nonlinearities; 
however, just for argument's sake, we restrict ourselves to the attractive 
case pT] . 

The hysteresis effect is strictly close to the appearance of spontaneous 
symmetry breaking phenomenon (SSBP) related to saddle point nodes. 
In fact, for BECs with (r+l)-body interaction governed by equation ([2]) 
has been recently seen [22] that SSBT appears when the nonhnearity 
power r is bigger than rthreshoid- It is worth mentioning the fact SSBP 
is a rather important effect that arises in a wide range of physical 
systems modeled by nonlinear equations [23] . We would also mention 
the fact that hysteresis effects associated to bifurcations of stationary 
solutions are theoretically discussed for BECs in optical lattice under 
the effect of a Stark-like external field [24i] 

The n-dimensional linear Schrodinger equation with a symmetric 
double well potential has stationary states of a definite even (^9+ and 
odd-parity with associate nondegenerate eigenvalues A_|_ < A_. 
However, the introduction of a nonlinear term, which usually models 
in quantum mechanics an interacting many-particle system, may give 
rise to asymmetrical states related to SSBP. 

In the semiclassical limit has been proved that the symmetric stable 
stationary state bifurcates when the adimensional nonlinear parameter 
7] takes absolute value equal to the critical value 



(3) 




(4) 



r]* = T jr . 
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The parameter t] is associated with the couphng factor of the nonhnear 
perturbation by 

(5) r] ■=r]{t) = cgr+i{t)/uj 

and it is the effective nonhnear couphng factor, where co is the (half of 
the) sphtting between the two levels 

(6) ,, = i(A_-A+) 
and c is the constant given by 

Here, (ffi and (fi are the normalized right and left hand-side vectors 

LfR = (v3+ + V5-)/V2 

and 

usually named single-well states because they are localized on only 
one well. In fact, in the semiclassical limit (or also for large distance 
between the two wells) the splitting u is exponentially small, as h goes 
to zero, and the supports of the two vectors ipR and ipi don't overlap 
up to an exponentially small term.. 

By adopting the two level approximation then the wave function 
ilj{x,t) is a linear combination of the right and left hand-side vectors 

^ix,t) = aR{t)LpR{x) + aL{t)<fL{x) 

where we set 

ttR = pe^", Ql = qe'^, + = 1 . 

Defining the relative phase difference 6 = a — P and the imbalance 
function z = — g^, and rescaling the time as r = ut/h (hence, 
the linear beating period takes the value vr), then equation ([2]) can be 
written in the Hamiltonian form 

86 _m dz _ an 

Or dz dr dO 

with Hamiltonian function 

(1 + zr+^ + (1 - zY+^ 



V- 



2'-(r + 1) 

The energy functional S associated to the nonlinear Schrodinger equa- 
tion ([2]) and written in the two level approximation takes the form 
£ = Q — ^uTi, where fl = |(A_ + A4.) is the mean value between the 
two energy levels. 



HYSTERESIS EFFECTS IN BOSE-EINSTEIN CONDENSATES 



5 



We consider at first the case r < rthreshoid- Since rj takes negative 
values then the nonhnear ground state is a stable symmetric state for 
any \ri\ < r]* . At jr^l = 77* it bifurcates and we observe also an exchange 
of the stability properties: for |?7| larger than 77* then the symmetric 
stationary state becomes unstable and the new asymmetrical states are 
stable (see Fig. [T]-a). 

On the other side, for r > rthreshoid, then a couple of saddle-node 
bifurcations, associated to new asymmetrical stationary states, sharply 
appears when jr/l is equal to a given value 77"*" such that 77"*" < 77* [25j; 
then, for increasing values of |7/|, the two unstable solutions disappear 
at 1 7/ 1 =77* showing a subcritical pitch-fork bifurcation (see Fig. [2]-a). 

In such a scenario, that is the sharp appearance of new asymmetrical 
stationary solutions fully localized on a single well when r > rthreshoid-, 
a new relevant effect occurs: namely we expect to observe hysteresis 
effect when we adiabatically changes the effective coupling factor r] such 
that its absolute value moves from values less than 77* to values bigger 
than 77* and then it goes back to its initial value. To this end we 
consider a state that, in the (z, 6') -represent at ion, is initially close to 
the symmetric stationary state: that is its initial condition corresponds 
to Zo ~ and 6q ^ 0. 

In the case r < rthreshoid then, as \r]\ increases, the state remains 
close to the symmetric stationary state for any I77I < 7/*, at I77I = 7/* it 
makes experience of a bifurcation and it follows one of two branches 
for |?7| > 7/*. When |7/| decreases from values bigger than 7/* to values 
less than 77* then such a path is reversed and the state returns close 
to the initial symmetric stationary state when \ri\ returns to its initial 
value such that \ri\ < 7/* (see the path indicated by the arrows in Fig. 
[U-a). 

On the other side, if r > rthreshoid then the state still remains close 
to the initial stable stationary state {zq,9q) for any \ri\ < rj*, but at 
|?7| = 77* we don't have a smooth bifurcation and for \r]\ > rj* the state 
starts to oscillate around the stable asymmetric stationary solution lo- 
calized on only one of the two wells. As I77I decreases from values bigger 
than 7]* to values between r/"*" and rj* the previous path is not reversed. 
In fact, the state continues to oscillate around the stable asymmet- 
ric stationary solution until \ri\ reaches the value r]~^. Then, while rj 
is returning to its initial value it takes the value \ri\ = r]~^, for which 
the asymmetrical stable stationary states disappear, and the wavefunc- 
tion starts to exhibit a wide oscillating motion around the symmetrical 
stationary solution corresponding to 2; = 0. If we introduce a small 
damping effect then such oscillating motions are damped and the state 
will stay close to the symmetric stationary solution. In conclusion we 
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can see in such a scenario that an hysteresis effect appears for values of 
|?7| between ri~^ and rj* (see the path indicated by arrows in Fig. |2]-a). 

In fact, such an hysteresis effect becomes more evident by adding a 
damping term which forces the state to collapse to the ground state. 
Actually, in physical systems we should expect to take into account a 
certain amount of damping due to the incoherent exchange of normal 
atoms. In particular, an accepted model for damped two-components 
BECs has been introduced by [25] and it reads as 

where > is the damping constant, and 



de &H 2ZC0S9 7] , xr- /i 

— [(1 + 2;) - (1-2:, 



In Fig. [T]-b (for r = 1) and Fig. [2]-b (for r = 5) we plot the numerical 
solutions of this dynamical system, where we assume the initial condi- 
tion ^0 = and zq = 0.01 closed to the symmetric stationary solution, 
and for times r G [0, T] where T = 4000. The damping factor is chosen 
to be z/ = 0.5 and the time dependent function t] has the following form 

r/(r) = -1- 2 [l-|2r/r-l|], if r = 1, 

and 

ri{r) = -3 - 5 [1 - \2t/T - 1|] , if r = 5. 

As predicted by means of the previous analysis on the bifurcation of 
the stationary solutions, it appears that for r = 1 no hysteresis effect 
occurs; while, for r = 5 the hysteresis effect occurs for \r]\ between r]^ 
and 7]*. Oscillations of the state that occur when \ri\ becomes less than 
7]* are damped because of the damping factor. It is worth mentioning 
also the fact that the delay observed in the case r = 1, when the 
absolute value of r] becomes larger than the branch point rj*, is not 
a consequence of some hidden physical effects but it comes from the 
singularity associated to the branch point. In fact, for larger values 
of T then rj is almost constant around the branch point and this delay 
disappear. 



In summary, we have shown that in a damped BEG with a double- 
well trapping potential the nonlinear term coming from a (r + l)-body 
interaction, for r bigger than 4, may give rise to an hysteresis effect 
when the corresponding coupling factor adiabatically changes. The 
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Figure 1. In this figure we consider the case r = 1 < 
^threshold- In panel (a) we plot the bifurcation diagram 
where full lines represents stable stationary solutions and 
dot lines represent unstable states. The arrows repre- 
sent the "path" of a solution initially close to the sym- 
metric ground state. A state initially close to the sta- 
ble symmetric stationary solution makes experience of 
a bifurcation effect at jr^l = 77* = 2 and it follows the 
new asymmetric stable stationary solution on one of the 
branches for |r/| > r/*. When |?7| moves from values big- 
ger than ?7* to values less than rj* then the state returns 
to be close to the stationary symmetric state without ex- 
hibiting hysteresis phenomenon. In panel (b) we show 
the numerical solution of the equations ([7]) and ([8]) for 
r = 1, 1/ = 0.5, zo = 0.01 and ^0 = 0. 



modulation of the coupling factor can be performed by tuning the con- 
densate population by means of suitable external fields. Such an hys- 
teresis effect has not theoretically predicted by means of the well-known 
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Figure 2. In this figure we consider the case r = 5 > 
rthreshoid- A statc initially close to the stable symmetric 
stationary solution jumps to the new asymmetric stable 
stationary solution on one of the branches for \ri\ > t]* = 
6.4. When |?7| moves from values bigger than 77* to values 
less than rj* the state does not return to be close to the 
stationary state, but it exhibits an hysteresis phenomena 
for |?7| between rj* and ri~^. Broken arrow are associated 
to the damping effect which forces the state to collapse 
on the ground state after some damped oscillations. In 
the right panel we show the numerical solution of the 
equation ([7]) and ([8]) for r = 5, = 0.5, zq = 0.01 and 
^0 = 0. 

Gross-Pitaevskii equation, in which only binary contact potentials are 
considered. 
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